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The high- Tc glass model can be combined with the repulsive tt '-Hubbard model as microscopic description of the 
striped domains found in the high-Tc materials. In this picture the finite Hubbard clusters are the origin of the 
d-wave pairing. In this paper we show, that the glass model can also explain the critical currents usually observed 
in the high-Tc materials. We use two different approaches to calculate the critical current densities of the high-Tc 
glass model. Both lead to a strongly anisotropic critical current. Finally we give an explanation, why we expect 
nonetheless a nearly perfect isotropic critical current in the high-Tc superconductors. 
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I. INTRODUCTION 

The high-Tc glass model was introduced in 1987 jl|,|§] to describe superconductivity in the high-Tc cuprates (HTSC) 
. Originally the glass model was designed for s-wave symmetry of the superconducting wave functions . Whereas 
it was first argued ||6| , that the glass model is not applicable for d-wave symmetry, the experimental result of striped 
domains in the superconducting CuO-planes inside the high-Tc materials only gives rise to weak disorder. Thus 
the high-Tc glass model is applicable in this situation, too |@]. It was demonstrated [Q, that the high-Tc glass model 
including the tt '-Hubbard model as a microscopic description of the striped superconducting domains is able to 
explain e.g. the d-wave symmetry of the superconducting phase and the pseudogap above Tc in the density of 
states (DOS) [0,0. In the combined high-Tc glass and striped Hubbard model picture the stripes in the HTSC 
occur at least at the same temperature, at which the pseudogap in the DOS opens, because otherwise the striped 
Hubbard clusters, which are responsable for these gaps, do not exist. 

In this paper we will show, that within the high-Tc glass model the maximum critical currents found in high-Tc 
materials and their almost perfect isotropy in a-b-direction can be understood. Following two independent paths we 
calculate in the next section an upper bound for the critical current density jc'- first with a direct calculation from 
the high-Tc glass model [0 and second considering the extended Bean model Afterwards we offer an intuitive 
picture for the (observed) isotropy of the critical currents and explain, why the values of jc usually measured in high-Tc 
materials are smaller than these upper bounds. From these simple ideas follows a fabrication procedure, which could 
lead to a possible increase of jc in the high-Tc materials. 

II. CRITICAL CURRENTS IN THE HIGH-Tc GLASS MODEL 

In the high-Tc glass model a single superconducting CuO-plane is described as an array of striped domains, figure [^, 
[01 . The typical dimension of a single striped domain with constant superconducting phase inside the high-Tc glass 



model can be roughly estimated following e.g. Tsuei and Doderer ||5| and Tranquada et. al. [141 . We take for our 
calculation the values a w 100 A, 6 w 10 A, (e.g. 6 « 16 A in LaSrCuO jlll or « 23 A in YBaCuO H) and z « lOA 
with the directions of figure |^. This choice of a, 6, and z is not crucial for our considerations, where only their order 
of magnitude is important for the conclusions. 

In the following calculation we show that the glass-model leads to critical current densities, which are in agreement 
with the experiments. The starting point is the Hamiltonian for the glass model 0,^]: 



n = -J^cos{(j)t - (pk - A,k) ■ (1) 
{i,k) 



The phase factors Ai^k are given by 



■3 



A,k = ^ I Adl (2) 



with TotA = H ■ z, where z is the unity vector in z-direction, and with Ai^k — <ii,k ■ H we get 
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aiM = ^ ^ [Vk - • (3) 

In equation (|l|) up to (|^) $o is the elementary flux quantum, (f)i are the phases of the superconducting wave functions, 
J is the coupling energy between two clusters, {i, k) is the sum over all nearest neighbors, and Xi, Xk^ Hi, and yk are 
the coordinates of the center of gravity of the domains i respectively k. Finally H is the external magnetic field and 
A the corresponding vector potential [yj . 

In the high-Tc glass model only weak, "correlated" disorder was chosen in the framework of the square lattice [0. 
We should note, that this weak disorder in the high-Tc glass model can be described and Hi K i ■ h with 

the lattice constants a in x-direction and h correspondingly in y-direction [00]. 

To obtain the critical current density ji^k for a connection between domain i and k we consider the Maxwell equation 

rot Hi^k^ ji,k ■ (4) 
Next we assume a constant increase of the internal magnetic field analogous to the assumptions of the Bean model 

rotgg- = . (5) 
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The free energy F is given as F = — fc^TlnZ, where the partition function Z = X^j*} ^^P with the 

inverse temperature f3 = l/(fcsr) and the sum {$} is over all possible configurations of the phases <& — {(pi, . . . , (j>L'^) 
for a lattice with L x L sites {ks is the Boltzmann constant). The magnetization M of a sample is given by 
M = —1/V ■ dF/dH, where V is the volume of the sample. For the external magnetic field H = we obtain: 

M = T^^J {yk - y,) (sin (0, - (j)k)) (6) 

where (. . .) dennotes the thermal expectation value. 

The magnetization can be expressed in terms of the internal magnetic fields Hi^k 

Inserting equation (^) into equation (0), using the abbreviation / = J27r/$o and the volume V — La-Lb-z (geometries 
of a striped domain, figure 0) it follows 

^ ^2/ ^' ^^'^ {yk - Vi) (sin {cpi - (j)k)) ■ (8) 
Using the equations (H), (||), and (||) we get for the critical current density of the domains i and k: 

ji,k = — (sin {(pi - (pk)) ■ (9) 
az 

Here we made also use of yk — yi ~ 6, if i and k are nearest neighbors (n.n.) in y-direction and yk — yi ~ in the 
other cases. 

Now we average over all 2L^ bonds between n.n. in the lattice and introduce jo = I /az, where az is the area "used" 
by a single junction between two domains. We obtain: 

j = 7^^y^J^,k =307;^^ {sin {(p^-(pk)) ■ (10) 



The average of the sinus: 



2X2 •^^22.2 

{i,k) {i,k) 



(sin (</>, - (Pk)} - ^ E (<^' - -^fe)) (11) 

{t,k) 



can be obtained from the simulation (index sim) with 
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(i,fe) 

We now make following approximations: yk,sim — Vi.sim = 1, if « and fc are n.n. in y-direction, yk,sim — Ui.sim = 0, if i 
and k are n.n. in x-direction, and Xi^sim + Xk,sim ~ 2i neglecting the random placement of the domains in the glass 
model |]. 

With i = 0, 1,...,-L — 1 this leads to the magnetization 



With Msim — O.SAMsim from figure |^ it follows 
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(sin (0, - = -AM,,rn ■ (14) 

Ij — 1 

Here also as in the Bean model j — > jc leads to the critical state and the critical current density is given by: 

jc=joj^AMs^m (15) 

with jo = 27rJ/(az$o)- 

Now we calculate the critical current density in a second way. Applying Bean's formula jl^,|l^ in the anisotropic 
case jl^ is justified as we have periodic boundary conditions in y-direction leading to a very long sample in y-direction. 
Concerning the internal magnetic fields in the critical state we have the situation illustrated in figure ^. Therefore we 
have in the simulation in x-direction a sample size of double length. Thus the experimentally found magnetization 
Mexp is the magnetization corresponding to the triangle of figure ]^ 

With the Bean formula of the anisotropic case |l^ : 

Me.p = ^ (16) 

and therefore with the length of the sample I = 2{L — l)a (figure ^ we have 

Now we calculate the value of the "experimental" magnetization M^xp from Msim'- 

Mexp = y-^J 2 ^^'^ ~ ~ ^^>l 



with Xi ^ a - Xi^sim {xi^sim = 0, 1, L - 1) and yt ^ b ■ yi^sim iyi,sim = 0, 1, i - 1) we have 

1 27r 
L'^abz $0 



{i,k) 



Therefore with / from equation (||) and Mgim from equation ( [l^ ) and with AMgim = 2Msim (figure ||) equation ([l9| ) 
leads to 

Mexp = -/AM,,™ , (20) 

z 

which leads (with equation (0)) to the critical current density 

Jc^—-f^AMs^m ■ (21) 
az L — 1 



Using the abbreviation jo = I /{az) and measuring jc in A/cm^ jig ] we have finally 



joj^AMs^m (22) 



3 



as in our first calculation (equation (|l5|)) of the critical current density using Bean's assumption (equation (|^)) in the 
high-Tc glass model directly. 

Using the numerical values AMs^m « 6 and L — 1 = 15 from figure|^ and the experimental values a = 100 A — 10"^ m 
and z = loA = 10"^ m for the stripes in the HTSC [|5|j|] we get / « 4.2 • 10~''A. And therefore we have for the 
critical current density 

J, = — AM,,„, «1.7-10^^ , (23) 
az cm^ 

which is surprisingly close to the experimental values jc ~ 5- lO^A/cm^ at 4K and zero field for the best films and 
higher than jc in wires {jc < 10^ Acm~^, e. g. [§^,^). We want to note, that in figure ^ the AMgim was measured at 
T = 0.2 J which corresponds to 20 K for = 100 K. 

To calculate critical current densities with equation ([l5| ) resp. (|2^) we made use of Bean's assumption of a constant 
increase of the internal magnetic fields It is clear, that this is a rather crude approximation for the high-Tc glass 

model. Furthermore for the calculation of jc we used the approximation of the averaged sinus (equation (lll|)). This 
gives perhaps a wrong estimate of the true critical current density of the high-Tc glass model. 

To support our line of reasoning we turn back to the glass model and estimate jc in the glass model directly following 
Ebner and Stroud . From the knowledge of the size of the domains it is possible to obtain the critical current Jc 
of the glass model, which can be determined directly ||2o| , El| . |l2| 

Jc^^J (24) 
n 

with J = ksTc and Tc = 100 K. (fc^ is the Boltzmann constant, e the elementary electron charge and h the Planck 
constant.) Equation (M) gives an upper bound for the maximum current through a single Josephson junction in the 
high-Tc glass model [^0[, which is used to estimate the coupling energy J in the high-Tc glass model . To determine 
the critical current density jc in a-direction we consider in a single domain (in figure |l|) the area A in b-z-direction, 
through which the current flows, figure With A k, h ■ z Ki 10^^* m^ wc obtain the critical current density 

^•c = :^ = ^M£ = 0.667. 10«^ . (25) 
A n bz cm^ 

This is indeed an upper bound to jc in equation (|23| ) . Thus both calculations (equation (^|) and (0) ) of jc give the 
same formula for jc, which is lower than the upper boimd for the high-Tc glass model in equation (P5|). 



III. (AN)ISOTROPY OF THE CRITICAL CURRENT 

Next we consider the critical currents in the b-direction instead of the a-direction. First we determine the magne- 
tization (with / = (27r/$o)-'^) 

M = -AM,,™ « 260— , (26) 
z cm 

which is in good agreement with the experiments jl^, too. Note, that M is independent of the size of the domains 
in the a-b-plane. Repeating the above calculations we have for the critical current density in b-direction: 

]c = ~JaJ^AMs^m . (27) 

This is analogous to equation ([isl ) with a different jo — I /{bz), in which a is replaced by b. Inserting the experimental 
values b « a/10 = 10 A and z « 10 A we obtain 

jt«10.j: = 1.7.10«^ . (28) 
cm^ 

In b-direction the upper bound analogous to equation (|25| ) also exists. Of course this upper bound is now also about 
ten times larger than in a-direction. 

Thus the consideration of stripes in the high-Tc glass model leads to a strong anisotropy of about the factor ten for 
the critical current densities between a- and b-direction depending on the ratio a/b of the striped domains. But this 
relatively large factor was never reported in the experimental literature. 
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We want to note, that on the one hand the size of the stripes enters the calculation of jc reciprocally, but on 
the other hand the number of domains in one spatial direction influences jc, too. Thus a better knowledge of the 
size of the domains is desirable. Additionally a finite size scaling (or the simulation of different system sizes L) for 
the magnetization AMsim of the high-Tc glass model is necessary to calculate the critical current densities more 
accurately. 

Now we have two features (anisotropy of jc in a- and b-dircction and a jc, which is in the order of magnitude of 
the largest experimentally found jc in HTSC) of the high-Tc glass model, which do not agree with the experiment. In 
our opinion there are two possible mechanisms, which can lead to a small or vanishing anisotropy. The first one is an 
anisotropic coupling constant J (Ja ^ Jh) in the high-Tc glass model, which may be different in a- and b-direction. 
But this probably only reduces the anisotropy. And it is unlikely, that these anisotropic coupling constants will lead 
to an isotropic critical current density jc- 

We postulate therefore the existence of Weiss- type domains in the planes with the stripes, which are dominantly 
either pointing in a- or b-direction, figure The existence of the Weiss-domains on the other hand explains the 
relatively lower critical current densities jc found in experiments for single crystals and thin films and in particular 
their differences. 

The resulting weak links between the Weiss-domains (figure which have to be assumed to be (much) "weaker", 
than the weak links of the high-Tc glass model, restrict jc to lower values. Thus equation (|2^ ) and ( psi) are only upper 
bounds of jc, too. 

This possibility was first brought to our attention in private discussions with Tsuei and Doderer. In the light of the 
above estimated factor of a/5 « 10 this picture is in our view one possible explanation for the experimental situation. 
Crystals or thin films with a predominant direction of the stripes are therefore predicted to show this anisotropy. 
But the main problem are the weak links between the Weiss-domains. These obviously govern the final experimental 
measurement of the critical current density. This picture could also explain the differences from sample to sample. 
While Tc is for all samples (nearly) identical, jc is quite different. But in our theoretical calculations for the high-Tc 
glass model both Tc and jc only depend on the coupling J. 

Considering this picture it is clear, that improvements of jc rise in general the quality of the samples. Especially 
in the light of the possible fabrication of electronic devices, better high-Tc samples have to be obtained by removing 
the Weiss-domains and the weak links between these domains. Therefore it would be extremely fruitful to avoid the 
Weiss-domains or to restrict the influence of the resulting weak links. In this spirit we propose a receipt following 
from the fabrication of magnets. The magnets are cooled down in a strong magnetic field. For high-Tc-niaterials 
we propose a similar procedure. The cooling process should use the cooling schedule known from simulations. This 
schedule has been subject to extensive research in the field of (physical) optimization in particular the "simulated 
annealing" method [^ , ^ . 

Procedures developed in the optimization theory ^2 23 2j,|2^ may be directly transfered to the annealing of the 
high-Tc-materials. Furthermore the annealing has to be carried out in an electric field and/or while an electric current 
is flowing through the sample. 

We expect from these procedures a substantial increase of the quality of the high-Tc-materials. In particular the 
fabrication of electronic devises as e.g. transistors should benefit greatly from these ideas. In the fabrication of wires 
or polycrystals additionally the weak links between the superconducting grains limit the critical current densities and 
should be removed, too. 



IV. SUMMARY AND CONCLUSIONS 

In the high-Tc glass model the striped superconducting domains of the high-Tc materials are the domains 

of constant superconducting phases The size of the striped domains may be deduced by mainly theoretical 
considerations and magnetic measurements @]. These measurements and the corresponding theoretical framework 
were already known shortly after the discovery of the HTSC 0. Already there the size of the domains could be 
estimated correctly to be approximately 10'' 0|. This leads to the conclusion, that the glass behavior is not related 
to the early ceramic structure of the grains, that many but not all samples exhibit. 

From the high-Tc glass model the critical current density can be calculated following two different paths: first 
using the extended Bean model ||l6|,|l3| and hysteresis measurements of the high-Tc glass model and second directly 
from the definition of the glass model using Bean's assumption (ig) of the critical state model. Both approaches lead 
consistently to the same critical current densities jc and to a strong anisotropic critical current density for the a- resp. 
b-direction of jc, which is much larger than the anisotropy found experimentally. But they are close to the highest 
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measured eritical current densities in thin films. This strong anisotropy follows from the underlying striped shape of 
the domains. 

Taking the same type of approximation (following the Bean model) we obtain the same value of jc in the simulation 

as in the experiments. This is in our opinion a strong evidence in favor of the high- Tc glass model. 

Nonetheless it is puzzling, why this anisotropy in the a-b-plane was never reported for jc- We propose two extensions 
of the high-Tc glass model: first an anisotropic coupling 7^ Jb, which in our opinion can only reduce the anisotropy, 
and second the existence of Weiss- type domains, in which the predominant direction of the stripes is turned by 90°. 
The latter leads to nearly isotropic critical current densities and lowers jc on account of the " new" weak links between 
these areas with the same direction of the stripes. This explains, too, why Tc is nearly constant between difii'erent 
samples whereas jc is quite different. Removing these Weiss- type domains should lead to higher critical current 
densities jc, which may be done by simulated annealing of the HTSC (eventual in electric fields). 
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FIG. 1. Superconducting (SC) domains (stripes) in the CuO-planes described microscopically with the Hubbard model (HM). 
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FIG. 2. Geometry and size of a single striped domain in the high-Tc glass model. 
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FIG. 3. Hysteresis measurement of the high-Tc glass model, (analogous to figure 12 in [1]) 
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FIG. 4. Critical state in a single striped domain in the high-T — c glass model. 
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FIG. 5. Weak links between Weiss-domains. 
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